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ABSTRACT 

Using three-dimensional convection simulations it is shown that a sinusoidal variation of hor- 
izontal shear leads to a kinematic a effect with a similar sinusoidal variation. The effect exists 
even for weak stratification and arises owing to the inhomogeneity of turbulence and the pres- 
ence of impenetrable vertical boundaries. This system produces large-scale magnetic fields 
that also show a sinusoidal variation in the cross-stream direction. It is argued that earlier 
investigations overlooked these phenomena partly because of the use of horizontal averaging 
and also because measurements of a using an imposed field combined with long time averages 
give erroneous results. It is demonstrated that in such cases the actual horizontally averaged 
mean field becomes non-uniform. The turbulent magnetic diffusion term resulting from such 
non-uniform fields can then no longer be neglected and begins to balance the a effect. 
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1 INTRODUCTION 



Shear can play an important role in hydromagnetic dynamos. This 
is especially true of dynamos in astrophysical bodies that gen- 
erate magnetic fields on scales larger than the scale of the tur- 
bulent motions. Those types of dynamos are generally referred 
to as large-scale dynamos. Simul ations confirm that shear can be 
the sol e driver of dynamo action jBran denburg 2005; Yous ef et al.l 
[ l2008al lbl; lBrandenburg et al.ll2008ah . but there is no consensus as to 
what is the underlying mechanism for producing such large-scale 
fields. In addition to shear there are also other possible mecha- 
nisms producing larg e-scale magnetic field s. One important con- 
tender is the a effect (Steenbeck et al. 1966), which quantifies the 
effect of kinetic helicity on magnetic field gen eration. It can also 
be the sole driver of large-scale dynamo action dBrandenburd200 H ; 
iKapyla et al]2009bl) . 

When both shear and a effect act simultaneously, it be- 
comes even harder to identify the main drivers of large-scale dy- 
namo action. Although shear is generally believed to b e advanta- 
geous for large-scale dynamo action (e.g. [Tobias 2009]), it is con- 
ceivable that the two effects (a effect and shear) suppress each 
other at least partially. This is because, in the presence of strat- 
ificatio n or other inhomogeneities, shear itself can produce an a 
effect dRogachevskii & KleeorirJ |2003| ; iRadler & Stepanovl 120061 ; 
Kapyla et alj |2009al) . Its sign depends on the relative orientation 
of shear and stratification. The net a depends then on the pseudo 
scalar (2f2 + W) g, where 2fi and W are the vorticities associated 
with rotation and large-scale shear flow, respectively. 

The issue can be complicated even further if shear is not con- 
stant but has a sinusoidal profi le, for example jBrandenburg et all 
l200ll ; lHughes & Proctodl2009l) . Sinusoidal shear profiles are com- 



monly adopted in numerical simulations where all boundaries are 
strictly periodic. This has obvious computational advantages and is 
certainly easier to implement than the so-called shearing-periodic 
boundary conditions where cross-stream periodicity applies only 
to pos itions that follow the shear flow and are thus changing with 
time ( Wisdom & Tremaine 1988). In helical turbulence with shear 
there is the possibility of dynamo waves that propagate perpendic- 
ular to the plane of the shear. This is clearly borne out by simula- 
tions dKapvla & Brandenburdl2009r) . The propagation direction of 
the dynamo wave is proportional to the product HkW, where Hk 
is the kinetic helicity of the flow. When the shear is sinusoidal, the 
sign of W changes in space, so one obtains c ounter-propagating 
dynam o waves in the two halves of the domain ( Brandenbur g et al.l 
120011) . In the presence of helicity, there is also a turbulent pumping 
effect, whose effe ctive velocity is also in the direction of HkW 
jMitra et al]|2009l) . 

In the cases discussed above the turbulence is driven by a he- 
lical body force, which is clearly artificial, but it allows contact to 
be mad e with analyti c theories of dynamo action in homogeneous 
media (Moffatt 1978). A more realistic case is one where the tur- 
bulence is driven by natural convection in a slab with a temperature 
gradient in the vertical direction. Many of the features of dynamo 
action discussed above carry over to this case as well, but an addi- 
tional complication arises both from the fact that there are impene- 
trabl e walls and that the sign of kinetic helicity changes with depth 
(e.g. lBrandenburg et al.lll990l ; ICattaneo & HughesfeoOrjl) . 

In the present paper we deal with both aspects, but we fo- 
cus in particular on the effects of sinusoidal shear, where we ex- 
pect at least partial cancellation of the a effect when averaged 
over horizontal planes. We contrast our work with earlier results 
that used linear shear, implemented via the shearing-box approx- 
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imation jKapvla et all [2008). as well as the case with no shear 
jKapvla et alj|2009bf) . where only the a effect can operate. The 
conclusion from these studies is that in the simulation domain there 
is an a effect of the strengt h expected from kinematic mean-field 
theory jKapvla et al.ll2009al lbl). There is also a back-reaction of the 
magnetic field through the Lorentz force, and its strength varies de- 
pending on wheth er or not magnetic he licity is allowed to escape 
from the domain jKapvla et al J 120093) . Again, these aspects are 
now well understood using mean-field th eory. The new aspect here 
is the sinusoidal shear. In a recent paper, Hughe s & Proctod j2009h 
present results from convection simulations with rotation and large- 
scale shear and report the emergence of a large-scale magnetic 
field whose growth rat e is proport i onal t o the shear rate, similar to 
the earlier results of Kapyla et al. (2008). They also determine the 
q effect from their simulations using the so-called imposed-field 
method and find that a is small and unaffected by the presence 
of shear. From these results the authors conclude that the dynamo 
cannot be explained by a classical a 2 or aQ, dynamo. 

The interpretation of the results of iHughes & Proctod j2009h 
is potentially in conflict with that of IKapyla et al.l j2008h . In both 
cases, convection t ogether with shear w as found to produce large- 
scale fields, but in Kapyla et al. ( 2008) they ar e interpreted as be- 
ing the result of a conventional a effect while in lHughes & Proctod 
(2009) it is argued that they are d ue to another mechanism similar 
to the incoherent q-shear effect (vishniac & Brandenburg! 1 19971 ; 
ISokolovll99^ISilant'evl200d:|Proctorl2007l). or perhap s the shear- 
current effe ct (Rogac hevskii & Kleeorirj 12003. 2004). Moreover, 
Hughes & Proctor (2009) argue that the a effect is ruled out. 

At this point we cannot be sure that there is reall y a difference 
in inte rpret ations, because the system s considered by Kap yla et al.l 
j2008l) and lHugh es & Proctod j2009l) are different in at least two 
important aspects. Firstly, in Hughe s & Proctor! j2009l) there is no 
density stratification, and since a is supposed to be pro portional to 
the logarithmic density gradient (Steenbeck et al. 1966) the result- 
ing a may indeed vanish. However, due to the impenetrable vertical 
boundaries, the turbulence is inhomogeneou s so that V In it rmB ^ 
0, which can also lead to an a effect (e.g. Gie secke et al J 12005 ). 
Here, Urms is the rms velocity of the turbulence. Secondly, the 
shear profile changes sign in the horizontal direction. Together 
with the vertical inhomogeneity this also produces an a effect 
jRogachevskii & Kleeorirj|2003l : iRadler & Stepano\]|2006l) . but its 
contribution is not captured by horizontal averaging and it partially 
cancels the a effect from rotatio n. This should be a measur able ef- 
fect which was not quantified in Hug hes & Proctor! {2009). Doing 
this is one of the main motivations behind our present paper. 

There is yet another important issue relevant to determining a 
in a system where the m agnetic Reynolds num ber is large enough to 
result in dynamo action jHubbard et al . 2009). Obviously, any suc- 
cessful a effect should produce large-scale magnetic fields. Given 
enough time, this field should reach saturation. By employing a 
weak external field one might therefore measure a at a saturated 
level. Depending on boundary conditions, wh ich were unfortu- 
nately not specified in IHughes & Proctod j2009l) . the saturation can 
result in a catastrophically quenched a effect. Furthermore, here 
we show that even in the absence of a dynamo the electromotive 
force from long time averages refl ects not only a due to th e uni- 
form imposed field as assumed by IHughes & Proctod j2009h . but 
also picks up contributions from the additionally generated nonuni- 
form fields of comparable magnitude. These caveats in determin- 
ing a with an externally i mposed field were k nown for some time 
jOssendriiver et alll2002l ; IKapyla et aUl2006l), but they have only 
recently been examined in detail jHubbard et al.l 2003) an d were 



therefore not addressed by Hughes & Proctod j2009l) . This gives 

another motivation to our study. 

Here we use a similar simulation setup as iHughes & Proctod 

j2009h and derive the a effect with the imposed-field method. We 
show that the value of a determined by the method of resetting 
the magnetic field after regular tim e intervals yields a substa ntially 
higher value than that reported by IHughes & Proctod (2009). Fur- 
thermore, we show that for a sinusoidally varying shear, also the 
q effect will have a sinuso idal variation in the hori zontal direc- 
tion, hence explaining why IHughes & Proctod j2009l) did not see 
the contribution of shear in their horizontally averaged results. 



2 THE MODEL 

In an effort to compare with the study o f lHughes & Proctod J2009), 



we use a Cartesian domain with L x 



5d and L z 



d with 



< z < d, where d is the depth of the convectively unstable layer. 
We solve the usual set of hydromagnetic equations 



OA 

— = U x B - ?7/xo J, 



Dlnp 
Dt 



= -V-U, 



— = --Vp + g-20, x U+ -J x B 

Dt p p 

+ -V -2ispS + -(U -C/ <0) ) 
p r 



(1) 

(2) 

(3) 



£l = _P V .[/+±V.^VT + 2,S 2 + ^J 2 , (4) 
Dt p p p 

where D/Dt = d/dt + U ■ V is the advective time derivative, 
A is the magnetic vector potential, B — V x A the magnetic 
field, and J = p^ 1 V x B is the current density, po is the vacuum 
permeability, rj and v are the magnetic diffusivity and kinematic 
viscosity, respectively, K is the heat conductivity, p is the density, 
U is the velocity, g — —gz the gravitational acceleration, and 
ft = fio(0, 0, 1) the rotation vector. The fluid obeys an ideal gas 
law p — pe(7 — 1), where p and e are the pressure and internal 
energy, respectively, and 7 = cp/cy — 5/3 is the ratio of specific 
heats at constant pressure and volume, respectively. The specific 
internal energy per unit mass is related to the temperature via e = 
c\rT. The rate of strain tensor S is given by 

S., = + Uj,i) - \5ijV ■ U. (5) 

The last term of equation l|3j maintains a shear flow of the form 



U — Uq cos 



2tt(x — Xo) 



(6) 



where Uo is the amplitude of the shear flow, xo ~ —L x /2 is the 
position of the left-hand boundary of the domain, and r is a re- 
laxation time. Here we use a r = 20-y/ d/g which corresponds to 

roughly 3.5 convect ive turnover times. 

In their study, IHughes & Proctod j2009l) use the Boussinesq 
approximation and thus neglect density stratification. Here we use 
the PENCIL Codi|3 which is fully compressible . Howe ver, in order 
to stay close to the setup of IHughes & Proctoil |2009) we employ 
a weak stratification: the density difference between the top and 
the bottom of the domain is only ten per cent and the average Mach 
number is always less than 0. 1 . Hence the effects of compressibility 



http://pencil-code.googlecode.com 
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Table 1. Summary of the runs. Here Ma = U Tms / y dg, where C/ rms is 
the total rms velocity including the shear flow, Mao = «rms/ \fdg, and 

Brms = firms /fieq, where fi.q = 

Co as 2.3, and Ra = 10 5 in all runs. 



18. 



Run 


Ma 


Ma/Mao 


Sh 




Dynamo 


AO 


0.028 


1.00 


0.00 




no 


Al 


0.027 


0.98 


0.07 




no 


A2 


0.028 


1.01 


0.14 


0.70 


yes 


A3 


0.039 


1.42 


0.36 


1.15 


yes 


A4 


0.063 


2.28 


0.72 


1.97 


yes 


A5 


0.096 


3.47 


1.45 


3.99 


yes 



are small. The stratification in the associated hydrostatic initial state 
can be described by a polytrope with index m = 1. Unlike our 
previous studies (e.g. Kapy la et alj2008l ), no stably stratified layers 
are present. 

The horizontal boundaries are periodic. We keep the temper- 
ature fixed at the top and bottom boundaries. For the velocity we 
apply impenetrable, stress-free conditions according to 



d z U x = OMy = u, = o. 



For the magnetic field we use vertical field conditions 



that allow magnetic helicity to escape from the domain. 



(7) 



(8) 



2.1 Units, nondimensional quantities, and parameters 

Dimensionless quantities are obtained by setting 

d = g = po = c P = p = 1 , (9) 

where po is the density at 2 m = -^d. The units of length, time, 
velocity, density, specific entropy, and magnetic field are then 



(10) 



[x]=d, [t] = y/d/g , [U] = y/dg , 
[p] = Po , [s] = c P , [73] = y/ dg pajio 



The simulations are controlled by the following dimensionless pa- 
rameters: thermal and magnetic diffusion in comparison to viscos- 
ity are measured by the Prandtl numbers 



Pr 



Pm = 



Xo 



(11) 



where \o = K/(cppa) is the reference value of the thermal dif- 
fusion coefficient, measured in the middle of the layer, z m , in the 
non-convecting initial state. We use Pr = 0.6 and Pm = 2 in 
most models. Note that H ughes & Proctoil d2009h use Pr = 1 and 
Pm = 5, but based on earlier parameter studies dKapyla et alj 
l2009al lcT) we do not expect this difference to be significant. The 
efficiency of convection is measured by the Rayleigh number 



Ra: 



gd 4 / l_ds 

v\o V c P dz 



(12) 



again determined from the initial non-convecting state at z m . The 
entropy gradient can be presented in terms of logarithmic tempera- 
ture gradients 



1 ds s 
cp dz j 



Hp 



(13) 
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Figure 1. Root mean square value of the total magnetic field as a function 
of time for the runs listed in Table[T] 



with V = (91nT/ainp) Zm , V ad = 1 - I/7, and Hp being the 
pressure scale height at z = z m . 

The effects of viscosity and magnetic diffusion are quantified 
respectively by the fluid and magnetic Reynolds numbers 



Re 



Rm 



ijkf 



PmRe, 



(14) 



where u m a is the root-mean-square (rms) value of the velocity 
taken from a run where U =0, and kf = 2n /d is the wavenum- 
ber corresponding to the depth of the convectively unstable layer. 
The strengths of rotation and shear are measured by the Coriolis 
and shear numbers 

2fi „, S 



Co = 



Sh = 



fcf' 



(15) 



where S = 2-kUq/L x . 

The size of error bars is estimated by dividing the time series 
into three equally long parts. The largest deviation of the average 
for each of the three parts from that over the full time series is taken 
to represent the error. 



3 RESULTS 

3.1 Dynamo excitation 

We fi rst set out to reproduce the results of iHughes & ProctoJ 
(2009). To achieve this, we take a run with parameters close to 
theirs which does not act as a dynamo in the absence of shear 
(Sh = 0). For this baseline simulation we choose the parameters 
Rm ~ 18 and Co gj 2. 3. We then follow the same procedure as 
IHughes & Proctoil ( l2009h and gradually increase Sh whilst keeping 
all other parameters constant (Table Q3 and determine the growth 
rate A of the magnetic field. 

The time evolution of the rms-value of the total magnetic field 
from our set of runs is presented in Fig. Q] We find no dynamo 
for Sh = and for weak shear with Sh = 0.07, the growth rate 
of the field remains virtually the same as in the absence of shear. 
This can be understood as follows: imposing large-scale shear via 
a relaxation term effectively introduces a friction term for U y in 

places where U — U 7^ 0, hence lowering the Reynolds number 
somewhat. However, as the same relaxation time ru r msfcf ~ 3.5 is 
used in all runs with shear, we are confident that these runs can be 
compared with each other. As the shear is increased beyond Sh = 



© 0000 RAS, MNRAS 000. 000-000 



4 P. J. Kapyla et al. 




0.1 



1.0 




Sh 





Figure 2. Growth rate A of the total magnetic field, divided by the shear 
rate S as a function of Sh. 



0.07, the growth rate first increases roughly directly proportional to 
the shear rate S (Fig.O. However, for Sh > 0.72 the increase of 
the growth rate slows down similarly as in several previous studie s 
dYousef et al]|2008bl : iKapyla et al.ll2008l ; iHughes & Proctorl2 009 ) . 



B J B ~ 





x/d 



Figure 4. Magnetic field component B y averaged over the saturated state 
in time and over the {/-dimension from Runs A2-A5. 



3.2 Field structure 

In earlier studies where a homogeneous shear flow was used, the 
large-scale magnetic field in the saturated state was non-oscillating, 
showed little dependence on horizontal coordinates, and could 
hence be well represented by a horizontal average iKapyla et al.l 
2008). However, in the present case with sinusoidal shear, the field 
structure and tempora l behaviour can in princip le be more com- 
plicated. Furthermore, Hughes & Proctor] d2009h do not comment 
on the field structure in their study. In fact, the only evidence of a 
large-scale field in their paper is given in the form of spectra of the 
magnetic field. 

We find that in our simulations the large-scale field is non- 
oscillating. It turns out that the magnetic field shows an interest- 
ing spatial dependence. In Fig. [3] we show visualizations of the 
structure of the B y component from the runs with the weakest 
(Sh ~ 0.14) and the strongest (Sh w 1.45) shear in which dy- 
namo action was detected. In both cases it is clear that the strong 
large-scale fields are concentrated to one side of the computational 
domain whereas the other side of the box is almost devoid of strong 
coherent fields. This behaviour is even more striking when the field 
is averaged over y and t; see Fig. [4] In the next section we show 
that the region of strong large-scale fields coincides with the region 
where the a effect is strongest. 



3.3 a effect 

The origin of large-scale magnetic fields in helical turbulence is 
com monly attributed to the a effect in turbulent dynamo theory 



(e.g. 



Moffatil l 19781 ; iKrause & Radier1ll980l ; [Riidiger & Hollerbacb] 



12004 ). Result s for convection simu lations, making use of the test- 
field method (Kiip vla et alj |2009b). suggest that the a effect does 
indeed con tribute to large-scal e dynamo action in simulations pre- 
sented bv lKapyla etail d2008l) . However, it was also shown that, 
in order to fully explain the simulation results, add itional contri- 
butions from the shear-current and CI x J effects (Radler 1969) 
appear to be needed. 



On the other hand, IHughes & Proctorl (2009) claim that 
in their setup the a effect is small, unaffected by shear, and 
thus incapable of drivin g a large-scale dynamo. The setup of 
IHughes & Proctorj J2009T) is based on the Boussinesq approxima- 
tion whereby stratification is not present in their system. However, 
the impenetrable vertical boundaries also generate an inhomogene- 
ity, which, in a rotatin g system leads to an a effect of the form 
dSteenbecketalJ 19661) 



at ] = Qi(G ■ + a 2 (G i n j + OA), 



(16) 



where G; denotes the inhomogeneity and CI is the rotation vector. 
In Boussinesq convection with rotation the kinetic helicity and thus 
the a effect are antisymmetric around the midplane of the layer. 
In such cases it can be useful to average over one vertical half 
of the layer to obtain an estimate of a. We note that mean-field 
dynamo models have shown th at the details of the a profile can 
also play a significant r ole (e.g. iBarvshnikova & Shukurovlll987l : 
IStefani & Gerbethll2003l) . In what follows, we show in most cases 
the full profile of a and present average s over the upper half of th e 
domain only when comparing directly to Hughes & Proctor (2009). 
Since the simulations in the present paper are weakly stratified, 
only minor deviations from a perfectly symmetric profile can be 
expected to occur. 

Adding a shear flow of the form presented in equation ® 
produces large-scale vorticity W z oc sin J, where a: is a shifted 
and rescaled x coordinate with x = 27r(a; — xp)/L x . Such vor- 
ticity leads to an a effect (see, e.g. Rogache vskii & Kleeorinl2003l : 
iRadler & St epanov 20q1), 



„(W) 



qi(G • W)6ij + ot2{G t Wj + GjWi 



(17) 



which, in the present case, leads to a yy oc sin x . Thus, when both 
rotation and shear are present, a — a(x, z) is a function of both x 
and z. 

In order to measure the a effect, we impose a weak uniform 
magnetic field Boe y , with Bo ~ 4 • 10~ B eq , and measure the 
response of the relevant (y) component of the electromotive force. 
Our a is then obtained from 
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Sh-0,14 
t= 4000.0 




Figure 3. Magnetic field component B y in the saturated state from two runs with weak (left panel, Sh ss 0.14, tu rms fcf PS 700) and strong shear (right 
panel, Sh si 1.45, £u r msfcf ~ 350). The sides of the boxes show the field at the periphery of the domain whereas the bottom (top) panel depicts By from 
z = 0.05d (z = 0.95d). 
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Figure 6. The contribution of shear to the a effect according to equa- 
tion 12 It . averaged over the {/-coordinate and time for Runs A1-A5. 



x/d 



Figure 5. The coefficient a, averaged over the {/-direction and time for 
Runs A0-A5. 



a = a yy = £ y /B . (18) 

In contrast to the study o f lHughes & Proctodd2009h . we do not usu- 
ally allow the field that is generated in addition to Bo to saturate, 
but reset it after a time inter val At m 10 tu TUl s k{. Su ch a proce- 
dure was fir st introduced by Ossendri iver et al.l |2002) and it was 
used also in iKapvla et all d2006h to circumvent the complications 
that arise du e to the additionally g enerated fields. A more system- 
atic study of iHubbard et al.l d2009h showed that only if At is not 
too long, the kinematic value of a can be obtained if there is a suc- 



cessful large-scale dynamo prese nt in the system. Howeve r, in the 
present study and also in that of lHughes & ProctoJ (2009) there is 
no dynamo in the runs from which a is computed. We find that it 
is still necessary to use resetting to obtain the correct value of a 
even in the absence of a dynamo. However, we postpone detailed 
discussion of this issue to Section [l!4l 

Our results for a from runs with constant rotation and varying 
shear are shown in Fig. [5] We find that in the absence of shear, a 
is a function only of z and has a magnitude of about 0.6ao, where 
ao = l^irms is a reference value, and firms is taken from a ran with 
Sh = 0. When shear is introduced, a increases (decreases) in the 
regions of the domain where sinx > (sinx < 0). However, for 
strong shear, the contribution to a from shear no longer appears to 
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a/a c 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 




Left half (q>0) 



Right half (q<0) 

■4 4 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Sh 

Figure 7. The coefficient a averaged over the upper half (0.5d < z < d) 
of the domain from the left (x < 0, solid line) and right (x > 0, dashed 
line) sides of the box. 



be symmetric around x — 0. This can be understood in terms of the 
shear parameter 



q = 

where 



dU 



(0) 



d.r 



-M 



dU 



■(()) 



Ox 



= S sin x. 



(19) 



(20) 



The flow is linearly unstable for q > 2 (Rayleigh instability cri- 
terion). Although the maximum value of q in our simulations is 
about 1.25 it is clear that for Sh > 0.36 (with \q\> 0.31), the pro- 
file and the magnitude of a are no longer significantly affected by 
the increasing shear. In order to illustrate this we compute the con- 
tribution of a due to shear from runs with Sh 7^ by subtracting 
the a that was found in the absence of shear using 



(W) 



(21) 



where is the a obtained from Run AO with no shear but only 
rotation. The results are shown in Fig. [6] and clearly show that for 
small Sh ( < 0.14), the shear-induced a shows a sinusoidal varia- 
tion as a function of x. For larger shear the profile of oS w ' is no 
longer antisymmetric around x — 0. This could reflect the asym- 
metry of the results for q > (— L x /2 < x < 0) and q < 
(0 < x < L x /2), that was found earlier by Snellman et al. (2009) 
in a somewhat different context of forced turbulence under the in- 
fluence of rotation and shear. They found that the Reynolds stresses 
were significantly different in setups with different sign of Sh or q, 
and that this asymmetry became more pronounced when the mag- 
nitude of shear was increased. Similar behavior has been seen in 
the magnetohydrodynamic regime by Korpi et al. (2009) in the 
Reynolds and Maxwell stresses. 

We also observe that the magnitude of w- w ' does not signif- 
icantly change for Sh > 0.36. This could indicate that the a effect 
due to shear saturates and that a simple relation like equation d!7t 
is no longer valid. This is apparent from Figure [7] which shows 
q volume-averaged over the upper half of the domain separately 
for the left and right sides of the box. For weak shear (Sh<0.2) 
we find that a is linearly proportional to shear. For Sh<0.4 the 
values of a on both sides appear to saturate to constant values. 
The results thus imply that the coefficients ot\ and a 2 in equa- 



Table 2. Summary of runs with and without resetting with varying Bo. 
Run Bl corresponds to Run AO in Table[T] Co ?s 2.3, Sh = 0, and Ra = 
10 5 in all runs and the imposed field in normalised form is given by Bo = 

-Bo / -Bcq ■ 



Run 


Rm 


Bo 


a/ao 


Resetting 


Bl 


18 


4 ■ 10 


0.39 ± 0.05 


yes 


B2 


18 


0.04 


0.36 ± 0.03 


yes 


B3 


18 


0.11 


0.37 ±0.05 


yes 


B4 


18 


0.39 


0.63 ±0.21 


yes 


B5 


18 


1.25 


0.25 ±0.10 


yes 


B6 


18 


4.47 


0.06 ±0.05 


yes 


CI 


18 


4 ■ 10- 5 


0.09 ±0.06 


no 


C2 


18 


0.04 


0.12 ±0.09 


no 


C3 


18 


0.12 


0.09 ±0.02 


no 


C4 


18 


0.37 


0.12 ±0.05 


no 


C5 


18 


1.27 


0.08 ± 0.03 


no 


C6 


18 


2.22 


0.06 ±0.01 


no 


C7 


18 


4.10 


(1.10 ±0.34) ■ 10~ 3 


no 


Dl 


30 


4 ■ 10~ 5 


0.36 ±0.03 


yes 


D2 


30 


4 ■ 10~ 5 


-0.03 ±0.23 


no 



ti on (TP7l> should depend on W when shear is strong. We note that 
in Hughes & Proctor (2009) also larger values of shear were used. 
The large vortex seen in the velocity field in their Figure 3 indi- 
cates that some of their runs with strong shear could indeed be in 
the Rayleigh-unstable regime. 

With the present data we cannot ascribe the appearance of the 
large-scale dynamo solely to the a effect. However, the coincidence 
of regions of strong magnetic fields and large a suggest that the 
a effect is indeed an important ingredient in generating the large- 
scale fields. 



3.4 Importance of resetting 

It has previously been demonstrated that the imposed field method 
can yield misleading results if a successful large-scale dynamo 
is oper ating in the system and long time averages are em- 
ployed faubbard et ai]|2009l) . In this case, unexpectedly low val- 
ues of a could be explained by the fact that the system is 
already in a saturated state. However, many papers have re- 
ported smal l values of a also for systems that do not act as dy- 
namos (e.g. ICattaneo & Hughes! I200f3 ; iHughes & Cattanedl2008l: 
iHughes & Proctorl I2009TT These results in apparent contradic- 
tion w ith those of lOssendriiver et al. 1 42002b : iKapyla et al.l J2OO6I 
2009a) who use either the imposed field method with resetting or 
the test field method. In these cases the syst ems must be in a truly 
kinematic state. Thus, the explanation of iHubbard et alj ( 2009) 
does not apply. The purpose of this section is therefore to resolve 
this puzzle. 

We begin the investigation of this issue by performing two 
sets of simulations where we study the dependence of a, as mea- 
sured using equation J 1 8b . on Bo with runs where the field is being 
periodically reset or left to evolve unimpeded (Sets B and C, see 
Table|2j. We take Run AO with Rm ~ 18 and no shear as our base- 
line and vary Bo /B cq in the range 4 ■ 10~ 5 ... 4. Our results for a, 
defined as the volume average over the upper half of the box, 



2 

IT 



Bo 



dz, 



(22) 
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Figure 8. Coefficient a according to equation {22} as a function of Bq from 
runs where the field is either being reset (Set B, solid line) or left to evolve 
on its own (Set C, dashed line). The dotted lines show fits to a quenching 
formula given in the legend where we use the coefficients q\ = 0.4 (0.1) 
and <J2 = 0.5 (0.2) in the upper (lower) curve. The diamonds on the left of 
the vertical axis indicate the values of a for Bo/B C q fa 4 ■ 10 . 
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are shown in Fig. [8] We see that, with the exception of the strongest 
Bq case in Set C, the results for both sets are in accordance with a 
simple quenching formula 

°~ l + g 2 (Bo/B cq ) 2 ' 1 ] 

where gi and q-i are constants which we use as free parameters 
in the fitting. We find that the value of a for weak fields is con- 
sistently four times smaller in the cases where no resetting is per- 
formed. The values of a in the range Bo/B e q « 0.04 ... 1 are 
essentially the same as those made for our standard imposed field 
strength Bo/B cq ~ 4 • 10~ 5 (see also Table[2]l. This suggests that 
the values of a in this range represent the kinematic stage and that 
the factor of four between the results in the different sets arises 
from the additional inhomogeneous mean magnetic fields gener- 
ated in the cases where no resetting is performed. 

This is demonstrated in the uppermost panel of Fig. [9] where 
the additionally generated horizontal magnetic fields, averaged 
over time and horizontal directions, are shown from Run CI. The 
origin of these fields can be understood as follows: the imposed 
field Boe y induces a z-dependent electromotive force in the y di- 
rection, i.e. £ y (z). This leads to the generation of an x component 

of mean magnetic field via B x (z) = . . . — £ y , z which, on the other 
hand, induces a z dependent electromotive force £ x (z) and hence 

B y (z) — . . . + £ X ,Z' Since these additional fields are functions 
of z, mean currents J x {z) = —B VtZ and J y (z) = B XtZ are also 
present. We emphasize that these fields arise due to the presence of 
an imposed field and decay if the imposed field is removed. 

It is now clear that a cannot be determined using equation d 1 8fc 
in this situation because the electromotive force picks up contribu- 
tions from the generated fields according to 

£ y (z) = a(z)[By(z) + B ] - Vt {z)J y (z). (24) 

Here we omit the off-diagonal components of and rjijk whose 
influence on the final result is marginal. Since the magnetic fields 
are weak, a and r/t can be considered as the kinematic values. We 
use here a as determined from Run B 1 (imposed field with reset- 
ting) and rjt obtained from a corresponding test field simulation 
(see the middle panel of Fig. [9jl where the test fields have a sin kz 



Figure 9. Top panel: horizontally averaged horizontal components of the 
magnetic field from Run CI. Middle panel: vertical profiles of ct(z) from 
the imposed-field method (solid line) and test field calculation with k = ki 
(dash-dotted line), and rjt(z) (dashed line). Bottom panel: j/-component of 
the electromotive force (solid line) compared with ctB y — rjtJy (dashed 
line), and aB y (dash-dotted line). 



dependence on z with k/ki = 1 and fci = 2n/d. For more details 
ab out the test field meth od in the context of convection simulations 
see lKapvla et al .(2009a). We normalise the turbulent diffusion with 
a reference value rjto = litrmsfcf -1 - The bottom panel of Fig. [9] 
shows that equation ( 1241 ) with these z-dependent coefficients gives 
a good fit to the simulation data of £ y from Run CI when the actual 
mean magnetic fields are used. The diffusion term in equation d24b 
has a noticeable effect only near the boundaries where the current 
is also largest. These results demonstrate that the interpretation of 
the electromotive force in terms of equation < | 1 3b is insufficient if 
long time averages are used. 

A general comment is here in order. Near boundaries, as well 
as elsewhere in the domain where the scale of variation of the 
mean field becomes comparable with the scale of the turbulent ed- 
dies, a simple multiplication with turbulent transport coefficients 
becomes inaccurate and one needs to resort to a convolution with 
integral kernels. The kernels can be obtained via Fourier trans- 
formation using the test-fi eld results for different wavenumbers 
(Brandenbu rg et alj2008bl) . In the present paper we have only con- 
sidered the result for the wavenumber k = 2ir/L z . This is also the 
case for the r)t shown in the middle panel of Fig. [9] The a obtained 
from the test-field method has a more nearly sinusoidal shape, but 
with similar amplitude than the profile shown in Fig. [9] This con- 
firms the internal consistency of our result. 

Another facet of the issue is highlighted when the magnetic 
Reynolds number is increased from 18 to 30 (Runs Dl and D2, 
see Fig. |10t. The larger Rm value is very close to marginal for dy- 
namo action whereas the smaller value is clearly subcritical. We 
find that, if resetting is used, the kinematic value of a is indepen- 
dent of Rm in accordance with mean-field theory. The situation 
changes dramatically if we let the field evolve without resetting; 
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Figure 10. Time series of the coefficient a for Runs B 1 and D 1 (uppermost 
panel), CI (middle), and D2 (bottom). 



comes meaningless without the knowledge of turbulent diffusion. 
The situation has now changed through the widespread usage of 
the t est-field method to obtain values of rj t at the same time (see, 
e.g., iGressel et alj|2008h . Furthermore, we show that, if the mag- 
netic field is reset before the additionally generated fields become 
comparable to the imposed field, the kinematic value of a can be 
obtained by much shorter simulations and without the complica- 
tions related to gradients of B or statistic al convergence. These is- 
sues were a lready known for some time Jossendriive r et al1l2002l : 
lKapvlaetal.ll2006l) . but they have generally not been taken into 
consideration. 

Another new aspect is the si nusoidal shear that is expe cted to 
lead to a sinusoidal a profile (e.g.lRadler & Stepanov 2006). In the 
study o f Hue hes & Proctod (2009) a volume average of a over one 
vertical half of the domain is used, which averages out the contri- 
bution of a due to shear. We find that, in the absence of shear, a 
is approximately antisymmetric with respect to the midplane of the 
convectively unstable layer suggesting that the main contribution to 
a comes from the inhomogeneity due to the boundaries rather than 
due to density stratification. When sinusoidal shear is introduced 
into the system, an additional sinusoidal variation of a in the x di- 
rection is indeed present. When the shear is strong enough, the a 
profile is highly anisotropic. The maximum value of a is close to 
the exp ected one, ap = |M rmB , w hich is significantly higher than 
the a in lHughes & Proctod j2009h . 

We also note that the regions of strong large-scale magnetic 
fields coincide with the regions where the a effect is the strongest. 
This supports the idea that the a effect does indeed play a signifi- 
cant role in generating the large-scale field. 



see the two lower panels of Fig.QJJ] For Run CI with Rm ~ 18 we 
can still extract a statistically significant mean value of a although 
the scatter of the data is considerable. For Run D2 with Rm « 30 
the fluctuations of a increase even further so that a very long 
time average would be needed to obtain a statistically meaning- 
ful value. A similar convergence issu e has been encountered in th e 
studies bvlCattaneo & Hughes! j2006T) ; lHughes & Cattanej J2008T) ; 
Hugh es & Proctod d2009l) . However, as we have shown above, the 
interpretation of such values cannot be done without taking into ac- 
count the additionally generated fields and the effects of turbulent 
diffusion. 



4 CONCLUSIONS 

We use three-dimensional simulations of weakly stratified turbulent 
convection with sinusoidal shear to study dynamo action. The pa- 
rameters of the simulations are chosen so that in the absence of 
shear no dynamo is present. For weak shear the growth rate of 
the magnetic field is roughly proportional to the shear rate. This 
is in accordance with earlier studies. A large-scale magnetic field is 
found in all cases where a dynamo is excited. The strongest large- 
scale fields are concentrated in one half of the domain (x < 0), 
with a sign change close to x — and weaker field of opposite sign 

in the other half (x > 0) of the box. 

In an earlier study, [Hughes & P roctod J2009h investigated a 
similar system and came to the conclusion that the dynamo can- 
not be explained by aft or a 2 dynamos due to a low value of a 
determined using the imposed-field method. However, we demon- 
strate that their method where long time averages are used yields 
the kinematic value a only if additionally generated inhomoge- 
neous mean fields are taken into account. Hence, this analysis be- 
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